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Abstract 


Given  a set  of  n half-spaces  in  three-dimensional  space  , we  develop  an 

algorithm  for  finding  their  common  Intersection  in  time  O(nlogn).  The 

intersection,  if  nonempty,  is  presented  as  a convex  polyhedron.  The  algorithm 
• • 
is  suamarlzed  as  follows:  (1)  the  half-spaces  are  placed  in  two  sets 

depending  upon  whether  they  contain  or  do  not  contain  the  origin;  (li)  the 

half-spaces  in  each  of  these  sets  are  dualized  to  points,  and  the  convex 

hulls  of  the  dualized  sets  are  obtained  in  time  O(nlogn);  (Hi)  since  the 

half -space  intersection  is  nonempty  if  and  only  if  these  two  convex  hulls  are 

disjoint,  a separating  plane  is  found,  also  in  time  O(nlogn) ; (iv)  after 

applying  a linear  spatial  transformation  which  maps  the  separating  plane  to 

infinity,  the  convex  hull  of  the  union  of  the  two  transformed  convex  hulls 

is  the  transformed  intersection  of  the  half-spaces.  Since  the  latter  can  be 

found  in  time  0(n) , the  overall  running  time  of  the  procedure  is  O(nlogn).  A 

significant  consequence  of  this  result  is  that  a three-variable  linear,  or 

convex,  programming  problem  can  be  asymptotically  solved  faster  than  by  the 

Simplex  algorithm,  in  the  worst  case. 
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1.  Introduction 

To  find  the  conmon  Intersection  of  a set  of  n three-dimensional  half- 
spaces is  an  Important  problem  in  computational  geometry.  Perhaps  the  most 
familiar  application  of  this  problem  occurs  in  mathematical  optimization, 
specifically  in  three-variable  linear  and  convex  programming  [1].  In  fact,  as 
is  well-known,  a linear  programming  problem  in  the  variables  (x,y,z)  consists 
of  extremizing  a linear  objective  function  of  x,  y,  and  z subject  to  a set 
of  n linear  inequalities  in  these  variables.  These  inequalities  define  the 
set  of  the  feasible  solutions  to  the  problem.  It  is  well-known  that  this 
set  is  a convex  region  of  space,  which  is  the  intersection  of  the  half -spaces 
corresponding  to  the  inequalities;  it  is  also  known  that  the  extremum  solution 
occurs  at  one  vertex,  or  extreme  point,  of  the  region  of  the  feasible  solu- 
tions (feasible  region).  Thus  the  linear  programming  problem  can  be  solved 
by  determining  the  extreme  points  of  the  feasible  region,  and  evaluating  the 
objective  function  at  each  of  them.  The  same  considerations  entirely  apply 
to  convex  programming  with  linear  constraints  when  the  convexity  of  the 
objective  function  is  such  that  its  extreme  value  occurs  at  an  extreme 
point  of  the  feasible  region. 

The  most  successful  method  for  solving  these  optimization  problems  is 
the  Simplex  algorithm  [1].  With  three  variables,  the  Simplex  algorithm 
involves  traversing  a path  on  the  polyhedral  surface  which  bounds  the 
feasible  region,  until  the  extremizing  vertex  is  reached.  Since  the 
Simplex  algorithm  spends  time  0(n)  in  moving  from  one  vertex  to  the  next 
vertex,  and  in  the  worst  case  it  may  have  to  visit  0(n)  vertices  before 
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2 

terminating,  we  conclude  that  its  running  time  is  0(n  ). 

For  the  two-variable  version  of  these  problems,  M.  I.  Shamos  and 
Hoey  [2]  have  shown  that  the  intersection  of  n half -planes  can  be  computed 
in  time  O(nlogn),  thereby  obtaining  a method  for  solving  the  two-variable 
linear  programming  problem  which  is  faster  than  the  Simplex  algorithm  in 
the  worst-case.  Their  method  is  based  on  the  divide-and-conquer  technique, 
i.e.,  on  finding  in  time  0(n)  the  intersection  of  two  polygons,  which  are 
respectively  the  common  intersections  of  two  sets  of  approximately  n/2 
half-planes . 

Shamos  [3]  also  suggested  that,  if  one  finds  a fast  algorithm  for 

intersecting  two  polyhedra,  then  this  could  be  used  - by  the  divide-and- 

conquer  technique  - to  obtain  a fast  algorithm  for  finding  the  intersection 

of  a set  of  half-spaces.  In  turn,  the  latter  could  be  applied  to  the  linear 

programming  problem.  In  a companion  paper  [4],  we  have  described  an  algorithm 

for  intersecting  two  polyhedra,  whose  total  number  of  vertices  is  n, 

in  time  O(nlogn).  This  method,  if  used  as  a merge  technique  in  a divide- 

2 

and-conquer  approach,  would  yield  a running  time  O(n(logn)  ) for  finding  the 
common  Intersection  of  n half -spaces.  In  this  paper  we  show.  Instead, 
that  the  latter  problem  can  also  be  solved  in  time  O(nlogn) . Our  technique 
is  crucially  based  on  the  polyhedron  intersection  algorithm,  but  not  as  a 
merge  technique;  rather,  it  uses  it  to  transform  the  original  problem  into 
its  dual,  i.e.,  that  of  finding  the  convex  hull  of  a set  of  n three- 
dimensional  points,  which  is  known  to  be  solvable  (by  a divide-and-conquer 
technique)  in  time  O(nlogn) . 

The  computation  model  adopted  in  the  preceding  algorithms  was  a 
random-access  machine  using  real-number  arithmetic.  We  shall  also  assume 
this  model  in  the  future  discussion. 
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We  shall  now  precisely  formulate  our  problem:  Finding  the  intersection 
of  n three-dimensional  half -spaces  consists  of  finding  the  solutions  to  a 
set  of  n linear  Inequalities  of  the  form 

ailX  + ai2y  + ai3Z  + ai4  ^ 0 (*■  * l,...,n)  (1) 

where  x,  y,  and  z are  cartesian  coordinates  of  a solution  in  three-space 
and  where,  for  each  1,  the  a^,  a.^,  a^*  a^  are  real  numbers  which  are 
not  all  0. 

For  reasons  to  be  explained  later,  it  is  more  convenient  to  express  the 
points  in  homogeneous  coordinates  x^,  x^,  x^,  and  x^  so  that  x * x^/x^, 
y ■ and  z * x3^x4*  We  obtain  therefore,  the  related  set  of 

inequalities; 

ailxl  + ai2x2  + ai3X3  + ai4x4  * 0 (l  ’ * <2> 

from  which  we  can  obtain  the  solutions  to  the  original  set.  In  fact,  the 

solutions  to  (2)  form  a convex  set  which  may  be  separated  into  three  disjoint 

subsets  each  also  convex.  They  are  (i)  the  positive  points,  having  x^  > 0 

which  correspond  to  solutions  of  (1),  (ii)  the  equatorial  points,  having 

x^  * 0,  and  (iii)  the  negative  points  having  x^  < 0. 

The  set  of  solutions  to  (2)  will  be  described  by  a system  having  the 
following  form: 

1.  A minimum  subset  of  (2),  (i.e.,  all  redundant  inequalities  will 
be  deleted) , corresponding  to  the  faces  F of  a generalized  polyhedron 
d,  to  be  defined  later. 

2.  A minimum  set  V'  of  extreme  points  (x^^.x^x^) . The  set  of 
solutions  to  (2)  consists  of  the  linear  combinations  of  the  elements  of 
V1  with  nonnegative  coefficients.  The  three-dimensional  nonhomo genaous 
correspondents  of  the  elements  of  V'  form  the  set  V of  vertices  of  the 
generalized  polyhedron  d. 


3.  A data  structure,  called  a doubly-connected  edge-list  (see  [4]) 


incident  on  vertices  in  V and  the  cycles  of  edges  bordering  faces  in  F 


One  advantage  to  treating  solutions  to  the  set  (2)  of  Inequalities  is  that 


the  extreme  points  all  lie  within  the  range  of  the  variables,  while  the 


extreme  points  of  (1)  may  lie  at  Infinity.  A second  advantage  is  that 


as  we  shall  see,  the  duality  between  V and  F is  complete  in  the  case  of  (2) 


and  we  shall  be  able  to  take  full  advantage  of  it 


In  the  next  section  we  shall  briefly  discuss  some  geometric  notions 


which  are  instrumental  to  the  subsequent  developments 


As  is  well-known,  there  is  a convenient  interpretation  of  homogeneous 


coordinates  in  E if  one  views  (x. ,x,,x,,x, ) as  coordinates  in  four-dlmenslcaal 


space  E and  normalizes  them  by  multiplication  by  a positive  constant 


1,  i.e.,  to  points  on  the  surface  of  the 


unit  hypersphere  S with  center  at  the  origin.  The  positive  points 


lie  in  what  we  shall  call  the  positive  open  hemisphere  consisting 


of  points  in  S with  x,  > 0;  similarly  the  negative  points  lie 


in  the  negative  hemisphere  and  the  equatorial  points  lie  in  the  inter 


section  of  the  hypersphere  with  the  equatorial  hyperplane  x 


sitlve  set 


An  illustration  of  the  positive  and  negative  sets  in  two  dimensions 
The  arrows  indicate  the  half-planes  defined  by  (1) . 
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The  three-dimensional  space  containing  the  solutions  to  (1)  can  be 

4 

viewed  as  the  hyperplane  x^  * 1,  which  is  tangent  to  S , and  a projection 

from  the  origin  establishes  a one-to-one  correspondence  between  points  of 

3 

the  positive  hemisphere  and  points  of  E . Each  of  the  hyperplanes  of  (2) 

4 

passes  through  the  origin  of  E and  defines  a correspoidlng  half-space; 

4 

the  common  intersection  of  these  half-spaces  is  a (convex)  cone  C-  . The 

4 4 

intersection  of  C with  S is  a connected  domain  which  may  cross  the 

4 4 

equatorial  hyperplane;  in  the  latter  case,  points  of  C-  ("I  S in  the  positive 


hemisphere  will  be  a projected  to  points  in  the  positive  set,  whereas  the 

4 4 

points  of  U fl  s in  the  negative  hemisphere  will  be  projected  to  a set 

of  points  in  the  hyperplane  x^  * -1,  which  will  in  turn  map  to  the  negative 

4 

set  in  x^  * 1 by  a symmetry  with  respect  to  the  origin  of  E (see  figure 
1(b)  for  an  illustration  in  one  less  dimension).  This  combined  set  in  the 
hyperplane  x^  « 1 may,  thus,  consist  of  two  separate  unbounded  convex  sets, 
the  positive  and  negative,  as  illustrated  two-dimensionally  in  figure  1(a). 
We  call  this  the  hyperbolic  case.  If  all  points  of  the  solution  project 
into  the  positive  (or  the  negative)  open  hemisphere,  the  corresponding  set 

3 

in  E is  bounded  and  we  call  this  the  elliptic  case.  The  parabolic  case 

3 

occurs  when  there  is  a single,  but  unbounded  set  in  E , i.e.,  the  case 

4 

when  equatorial  as  well  as  positive  (or  negative)  points  occur  in  S . In 
all  three  cases  the  sets  of  points  form  what  we  shall  call  a generalized 
convex  polyhedron.  The  elliptic  case  alone  corresponds  to  a conventional 
polyhedron. 


6 


. 


Ia  the  sequel  we  shall  make  frequent  use  of  an  lnvolutory  transformation 

from  points  to  planes  and  vice-versa  called  duallzatlon.  Under  dualization 

the  coordinates  (§^»§2»53»§4.)  a P°*-nt  are  reinterpreted  as  the  coefficients 

defining  a half-space  §jX^  + §2X2  + £3*3  + £4*4  S Dualization  is  susceptible 

4 3 

of  two  intuitive  geometric  interpretations  - one  in  E and  one  in  E . In 
4 4 

E each  point  on  S is  viewed  as  the  terminus  of  a (unit)  vector  applied 

to  the  origin  and  its  dual  is  the  half-space  containing  this  vector  whose 

boundary  hyperplane  passes  through  the  origin  and  is  orthogonal  to  this 

vector.  Conversely,  each  half -space  through  the  origin  dualizes  into  the 

3 

corresponding  vector.  In  E a positive  point  at  distance  l from  the  origin 
is  mapped  into  a half-space  containing  the  origin  whose  boundary  plane  is 
at  distance  1 !l  from  the  origin  in  a direction  opposite  to  the  point.  If 
the  point  is  negative,  the  complementary  half -space  is  obtained  with  the  same 
boundary.  Conversely,  any  half-space  whose  boundary  plane  does  not  contain 
the  origin  has  a point  as  its  dual.  This  point  is  interpreted  as  being 
positive  or  negative  depending  upon  whether  the  half-space  contains  the 
origin. 

After  dualizing  the  system  of  inequalities  we  shall  need  to  apply  various 
algorithms  which  construct  the  convex  hull  of  a set  of  points,  and  follow  this 
operation  by  an  inverse  dualization.  Since  the  convex  hull  algorithms  in  their 
published  version  [5]  apply  to  the  conventional  or  elliptic  case,  we  shall 

r 

find  it  convenient  to  first  apply  an  invertible  linear  transformation  of 

4 3 

coordinates  to  (2)  in  E . This  transformation  will  have  the  effect  in  E 
of  moving  the  origin  to  a selected  point  and  hence  will  enable  us  to  reduce 


|i 
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the  dual  system  to  the  elliptic  case  as  desired.  The  transformation  will 

be  described  by  a nonsingular  4x4  postmultiplicative  matrix  T.  It  is 

desirable,  however,  in  order  to  avoid  unnecessary  coordinate  manipulations, 

4 4 

to  restrict  ourselves  to  linear  transformations  of  E which  map  S to  itself. 
These  transformations  are  readily  characterized.  Letting  the  row  vector  § 

4 (11 

denote  a point  on  S , we  know  that  £•£.'  =*  lv  , the  image  PT  of  £ must  be 

4 4 

on  S , whence  £T(£T)'  * £TT'j[/  ■ 1.  Since  £ is  arbitrary  on  S , the  latter 

holds  if  and  only  if  T*T'  * 1^,  the  4x4  identity  matrix,  i.e.,  T * ■ T*. 

This  is  the  characteristic  property  of  rotations : thus  we  shall  only  consider 
4 

rotations  of  E . It  is  worth  noticing  that  rotations  commute  with  dualization, 
that  is,  the  dual  of  the  image  (under  rotation)  of  a point  coincides  with 
the  image  of  the  dual  of  that  point. 


r 

I ^ 


B 

I 
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3 . Finding  the  Intersection 

We  shall  now  rewrite  (2)  in  matrix  form 

Ax'  i 0 (3) 

where  A * |ja^j  is  an  n X 4 matrix.  Let  A+,  Aq,  and  A be  the  matrices 
formed  from  the  rows  of  A for  which  a^  > 0,  a^  = 0,  and  a^  < 0 
respectively.  In  a later  discussion  we  shall  show  that, except  in 
degenerate  cases  which  can  be  reduced  to  a problem  of  lower  dimensionality, 
it  is  possible  to  find  a rotation  RQ  of  A which  eliminates  Aq.  Therefore, 
we  assume  that  in  (3)  such  a rotation  has  already  been  made  and  that  either 

ai4  > ® or  ai4  < ® f°r  each  *-• 

^ is  the  transpose  of  £. 
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Each  row  (ai I,ai2,a13,ai4^  A+  rePresents  a Plane  E-5  (or,  equivalently, 

4 

a hyperplane  through  the  origin  in  E ) . We  now  dualize  each  such  plane,  i.e., 

3 

we  construct  a point  (a^/a^.a^/8^*4^/*^)  in  E and  aPPty  the 

(D) 

algorithm  of  [5]  to  find  the  convex  hull  7+  of  these  points:  notice 

that  the  algorithm  is  directly  applicable  because  each  point  has  finite 
coordinates.  The  same  procedure  is  then  used  to  construct  a convex 
polyhedron  7^  from  the  matrix  A . We  may  now  apply  to  the  two  three- 
dimensional  polyhedra  7^  and  7^  an  algorithm,  described  in  [4],  which 
determines  whether  the  two  polyhedra  intersect  and,  if  so,  finds  the  inter- 
section, and  if  not,  finds  a separating  plane.  Actually  if  7^  fi  (7 ^ has  an 
interior  (i.e.,  nonzero  volume),  we  need  not  explicitly  construct  it,  since 
in  this  case  the  set  of  inequalities  (3)  is  inconsistent  and  has  no  solution. 

In  fact,  by  simple  reasoning  on  duality,  any  point  u ■ (ui,u2’u3^  in  the 
interior  of  7^  must  satisfy  + U2X2  + ^x^  + x^  > 0 when  x ■ (x^,x2,x.j,x^) 

is  any  solution  of  (3).  If  u is  also  in  the  interior  of  7^  then  it  must 
satisfy  -UjX^  - XJ2x2  ~ U3X3  - > 0 as  well,  so  we  conclude  that  (3)  has  no 

solution. 

If  7^  fl  7^  is  nonempty  but  has  no  interior,  the  problem  reduces  to 
one  of  lower  dimensionality  to  be  discussed  in  Section  4. 

Finally  suppose  that  7^  and  7^  do  not  intersect.  In  this  case  we 
find  a separating  plane  by  the  method  described  in  [4],  and  let 
p x + PjY  + P^2  + p^  ■ 0 be  this  plane,  with  the  signs  of  the  coefficients 
so  chosen  that  for  any  point  (a,b,c)  in  af)  we  have  p^a  + p2b  + p^c  + p^  >0. 


The  letter  "D"  as  an  apex  is  to  remind  us  that  we  are  dealing  with  the 
convex  hull  of  the  dual  set  of  points. 
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The  separating  plane  just  £ound  may  be  extended  to  a hyperplane 

4 

PjX^  + p2x2  + p^Xj  + p^x^  » o through  the  origin,  whose  normal  in  E is  the 
vector  (Pi»P2»p3»p4^ • We  now  wish  to  rotate  coordinates  in  E so  that  this 

3 

vector  points  toward  the  origin  of  E . In  other  words,  we  wish 

to  find  a rotation  R such  that  (Pj_,P2,P3,P4)R  - (0,0,0,K),  where  K is 

/2  2 2 2 

a positive  number  (actually  K ■ vT>2  + p2  + p3  + p4>  * A *ult»1»le  rotation 
is  easily  constructed  as  a product  RjR2R3»  where  each  causes  p^  to  go 
to  0.  For  example,  we  may  take 


f + 

' o 

1 -Pi^i  + P4 


0 0 

1 0 

0 1 

0 0 


The  elements  of  the  4x4  matrix  R can  be  computed  in  time  bounded  by  a 
constant  using  the  operations  of  arithmetic  including  square  root. 

We  claim  that  the  fourth  column  of  the  matrix  AR  consists  entirely  of 
positive  entries.  In  fact,  system  (3)  can  be  rewritten  as 

Ax'  - A(RR'1)x*  - ARR'x'  - AR(xR) ' * 0.  (3') 

Since  the  point  j - (P1,P2,P3,P4)  satisfies  atlpx  + ai2P2  + ai;Jp3  + at4p4  > 0 
for  each  1,  the  image  of  £ under  the  rotation, which  is  (0,0, 0,K), 
satisfies  a corresponding  set  of  strict  Inequalities  using  the  rows  of 
AR.  In  each  of  these  inequalities  the  only  nonzero  term  is  the  positive 
constant  K multiplied  by  the  fourth  entry  in  the  corresponding  row  of  AR, 
thereby  establishing  the  claim. 
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1 : 


There  is  no  need,  however,  to  construct  AR,  since  it  may  contain 
irrelevant  rows  which  we  want  to  suppress.  Instead,  for  each  vertex 

(ail/ai4’ai2/ai4,ai3/ai4)  °f  con8truct  a “w  (ail’ai2,ai3,al4)  of 

a new  matrix  A+m^  whose  rows  form  a subset  of  the  rows  of  A+.  Similarly, 

for  each  vertex  of  G ^ we  find  the  corresponding  row  in  A and  from  these 

rows  form  a new  matrix  A^ . 

Let  A*  - A^R  and  A*  - A^R;  the  set  of  inequalities 

[A*.  A*]'  5'  20  (4) 

is  equivalent  to  AR§'  2 o,  with  all  irrelevant  inequalities  deleted. 

3 

Since  (4)  is  satisfied  by  (0,0,0, 1),  the  origin  in  E , conceptually  the 

intersection  of  the  original  n half-spaces  can  be  found  by:  (i)  dualizing  the 

planes  corresponding  to  the  rows  of  [A*,A*],  (ii)  finding  the  convex  hull 

of  the  set  of  points  thus  obtained,  and  (iii)  dualizing  back  the  convex  hull. 

However,  this  sequence  of  operations  can  be  considerably  simplified.  The 

duals  of  the  planes  corresponding  to  the  rows  of  A^  and  A^  are  already 

available  as  the  vertices  of  the  polyhedra  G^  and  G ^ respectively.  Thus, 

to  obtain  the  analogous  polyhedra  G*^  and  G*^  for  A*  and  A*  respectively, 

we  rotate  the  polyhedra  G^  and  G^  (which  lie  in  the  x^  = 1 hyperplane) 

4 

by  R and  then  project  them  through  the  origin  of  E back  to  the  hyperplane 
x^  * 1.  In  practice,  this  means  multiplying  the  coordinates  of  each  of 
the  rotated  vertices  by  a number  which  makes  its  fourth  coordinate  1. 

The  two  polyhedra  G*^  and  G*^  are  disjoint  because  they  are 

4 

separated  by  a hyperplane  through  the  origin  of  E normal  to  the  R-transform 
of  the  vector  (0,0,0, 1).  We  therefore  may  construct  the  convex  hull  of 
the  union  G*^  U G*^  by  a single  application  of  the  "me*  a"  portion 


mu  l ai  llTl.1  1 1 ■n.lffaiw  m 
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of  the  algorithm  of  [5]  In  time  proportional  to  the  total  number  of 
vertices  In  the  two  polyhedra.  Let  us  call  this  resulting  polyhedron 

Now,  the  dual  3*  of  the  polyhedron  3*^  represents  the  set  of 
solutions  to  the  system  of  Inequalities 


where  the  coefficients  a*j  are  the  elements  of  the  matrix  A*  - AR,  and  we 
recall  that  all  the  coefficients  a*^  are  positive.  We  therefore  apply  the 
Inverse  rotation  R’*  to  the  vertices  and  faces  of  C7*.  Now  3*  lies  in 
the  hyperplane  » 1,  so  the  four-dimensional  rotation  R~*  Is  applied 
to  vertices  of  O of  the.  form  (v^,vj,v*, 1) . The  resulting  vertices 
(v^v^.v-j.v^)  must  be  renormalized  if  v^-  + 0,  and  this  Is  done  by 
multiplying  each  of  the  coordinates  by  the  positive  constant  1/ |v^ | . 

The  result  of  this  rotation  and  renormalization  is  a set  of  vertices 
and  associated  faces  which  we  can  represent  by  a doubly-connected  edge- 
list  and  which  we  shall  call  a generalized  polyhedron  3.  The  vertices  of 
this  generalized  polyhedron  with  fourth  coordinate  equal  to  1 are  the 
extreme  points  of  the  intersection  of  the  half-spaces. 

In  the  beginning  of  this  discussion  it  was  stated  chat  an  initial  rotation 
Rq  of  coordinates  could  be  found  which  would  make  each  entry  a^  in  the 
fourth  column  of  the  matrix  A different  from  0.  We  shall  now  describe  how 


12 


Rq  is  obtained.  As  before,  we  define  Aq  to  be  the  matrix  composed  of  rows 
from  A having  a^  ■ 0.  Inequalities  of  the  set  (1)  corresponding  to  these 
rows  appear  as 

*11*  + *12*  + *13z  2 0 • <-t"> 

3 

We  seek  to  find  a point  u -(u^.^.u^)  in  E which  strictly  satisfies  all 

these  inequalities.  Nbtice  that  the  set  of  Inequalities  (1")  defines  a 
3 3 

convex  cone  <3  in  E , whose  vertex  is  the  origin.  Except  in  degenerate 
3 

cases,  C-  intersects  either  one  of  planes  z - 1 and  z - -1,  since  it  cannot 

be  contained  in  the  region  -1  < z < +1.  Hence,  if  a point  u exists,  it  can 

be  found  as  follows.  Let  z * 1 in  (1")  and,  using  an  algorithm  due  to 

Shamos  and  Hoey  [2]  for  finding  in  time  O(nlogn)  the  Intersection  of  n 

half-planes,  find  the  polygon  of  solutions.  If  no  such  polygon  exists, 

then  set  z » -1  and  repeat  the  process.  In  one  case  or  the  other  there 

will  be  a polygon  P from  which  we  can  choose  _y  as  one  of  its  interior  points. 

4 

Now  we  augment  u » (ui»u2,U3^  to  form  a point  v * (u^u^u^M)  in  E , 
taking  M > 0.  Clearly,  v strictly  satisfies  all  the  inequalities  in 
AqV  ^ 0.  As  we  saw  in  the  derivation  of  R,  we  can  construct  a rotation 
Rq  such  that  vRq  ■ (0,0,0,Kq),  and  this  rotation  has  the  property  that  all 
entries  in  the  fourth  column  of  AqRq  are  positive. 

The  rotation  Rq  can  be  made  to  differ  very  little  from  the  identity 
transformation  by  taking  M very  large.  Thus,  we  take  M large  enough  so 
that  no  fourth  column  entries  in  A+  or  A change  sign  when  the  rotation  Rq 
is  performed.  Hence  all  entries  in  the  fourth  column  of  AIL  are  nonzero. 


I i 
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As  v«  noticed  construction  of  u can  be  done  In  time  O(nlogn) ; the 
construction  of  a rotation  R^,  can  be  done  In  constant  time,  as  previously 
discussed  In  a similar  situation. 

In  conclusion  we  will  replace  the  original  set  of  inequalities  (3)  by 


I 


0 

1. 

0 

n 

r. 


% 


(3”) 

To  summarize  the  preceding  discussion,  the  intersection  of  n half-speces, 
expressed  by  (3), can  be  obtained  as  follows. 

Step  1:  Find  a point  u in  the  interior  of  the  set  of  solutions  to  (l*) 
and  construct  a corresponding  rotation  Rq.  Replace  A by  ARq 
[This  can  be  done  in  time  O(nlogn)  using  the  half -plane 
intersection  algorithm  [2].] 

Step  2:  Find  the  convex  hull  0^  of  the  set  of  points  (*ii/*i4t»12^*i4,*i3^*i4^ 

for  all  i such  that  a..  >0.  Find  the  convex  hull  0 ^ of  the  set 

i4 

of  points  (ati/ai4»a12^ai4,ai3^ai4^  for  a11  1 such  that  ai4  < °* 

[This  can  be  done  in  time  O(nlogn)  by  using  the  convex  hull 

algorithm  described  in  [ 5 ] . ] 

Seep  3:  Using  the  algorithm  for  intersecting  polyhedra  [4],  test  if 

n 0 ^ has  an  interior.  If  so,  the  half-space  intersection 
is  empty;  if  not,  find  a separating  plane  pjX  + p2y  + p3z  + p^  - 0. 

[This  can  be  done  in  time  O(nlogn),  according  to  the  results 
reported  in  [4].] 

Step  4:  Find  a rotation  R of  so  that  the  separating  plane  is  mapped  to 

the  plane  at  infinity  [this  can  be  done  in  constant  time].  Find 
the  R-trans forms  of  the  polyhedra  and  and  project  them 

through  the  origin  back  to  the  - 1 hyperplane.  [This  requires 
time  0(n).]  Let  0*^  and  0*^  be  the  resulting  polyhedra. 


* 

- 


Step  S: 


Find  the  convex-hull  d*^  of  the  union  of  the  two  polyhedra 
and  d*^\  [This  is  done  in  tine  0(n)  using  the 
algorithm  of  [5] . ] 

Step  6:  Compute  the  R-1^1- transform  of  the  vertices  and  faces  of  the 
dual  <7*  of  d*^ . Then  renormalize  the  transformed  vertices 
to  obtain  d,  the  generalized  polyhedron  which  forms  the 
intersection  of  the  given  half-spaces,  [this  step  also 
requires  time  0(n).] 

In  conclusion,  we  see  that  the  preceding  algorithm  runs  in  time 
O(nlogn),  Steps  1,  2,  and  3 being  the  limiting  ones. 

To  eliminate  solutions  to  (2)  corresponding  to  negative  points,  one  can 
augment  the  set  (2)  with  the  Inequality  x^  * 0 before  applying  the  algorithm. 
The  remaining  positive  points  project  to  solutions  to  the  set  (1),  and  if 
there  are  equatorial  points  which  are  extreme  points  they  should  be 
retained  to  facilitate  the  description  of  the  set  of  solutions. 

4.  Degenerate  Cases 

One  degenerate  case,  which  may  arise  in  Step  3 of  the  main  algorithm, 
occurs  when  the  intersection  d^  D d ^ is  nonempty  but  has  no  interior, 
i.e.,  has  zero  voliane.  In  this  case,  the  algorithm  described  in  [4]  may  be 
used  to  find  a point  ^ <7<D)  n </»> . If  x - (xlfx2,x3,x4) 

is  any  solution  to  (3),  then  both  q^x^  + q2x2  + q^x^  + i 0 and 
-q^x^  - q2x2  - q3*3  • x^2  0 must  be  satisfied.  Hence,  the  set  of  solutions 
to  (3)  is  constrained  to  the  hyperplane  q^  + q2x2 '+  q3x3  + x^  • 0.  Thus, 

It  'it  it  It 

we  construct  a rotation  R such  that  (q^,q2,q3, 1)R  ■ (0,0, 0,K  ),  where  K 


is  a positive  constant 
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Now,  after  R is  applied,  the  set  of  inequalities  AR  x'  2 0 imply 

4t 

x^  ■ 0,  so  the  fourth  column  of  AR  is  superfluous.  Thus,  the  set  of 

Inequalities  AR*x'  ^ 0 may  be  replaced  by  a set  having  the  form  (1') 

whose  algorithmic  solution  has  already  been  described.  Hence,  to  obtain  the 

solution  to  (2)  one  need  only  rotate  the  result  by  (R*)"^. 

Another  degenerate  case  could  arise  in  Step  1 of  the  main  algorithm 

when  there  is  no  interior  point  u to  the  set  of  solutions  to  (1'),  i.e. 

to  AqX1  ^ 0.  In  this  case,  while  seeking  such  an  interior  point,  we  either 

discover  (a)  that  while  no  interior  point  exists  there  is  a point 

v * (v^.v^.Vj),  with  v3  * 1 or  v3  » -1  which  satisfies  (1'),  or  (b)  that 
3 

the  cone  <3  fails  to  intersect  either  plane  z • +1  or  i « -1. 

In  case  (a)  the  point  v must  lie  in  a plane  through  the  origin  which 

contains  all  solutions  to  (1*).  Suppose,  this  plane  has  the  equation 

* ★ * 
m^x  + m^y  + m^z  - 0.  We  then  choose  R so  that  (m^.n^.m^C^R  * (0,0, 0,K  ) 

for  some  positive  K*  and  solve  as  in  the  previous  case.  The  equation  of 

the  desired  plane  Is  obtained  as  a result  of  applying  the  algorithm  due  to 

Shamos  and  Hoey  [2], 

In  case  (b) , the  inequalities  (1')  imply  that  z * 0.  Hence,  the  third 
column  of  A is  superfluous.  The  set  of  inequalities  Ax'  ^ 0,  given  in  (3), 
may  be  replaced  by  a set  of  the  form  (1')  by  deleting  the  third  column. 

The  solution  to  (!')  is  obtained  as  before. 
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